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O ' Abstract 

in ' 

Q\ . We investigate cosmological structure formation seeded by topological de- 

fects which may form during a phase transition in the early universe. First we 
Q_i' derive a partially new, local and gauge invariant system of perturbation equa- 

tions to treat microwave background and dark matter fluctuations induced 
by topological defects or any other type of seeds. We then show that this 
system is well suited for numerical analysis of structure formation by apply- 
ing it to seeds induced by fluctuations of a global scalar field. Our numerical 
results cover a larger dynamical range than previous investigations and are 
complementary to them since we use substantially different methods. The 
resulting microwave background fluctuations are compatible with older simu- 
lations. We also obtain a scale invariant spectrum of fluctuations with about 
the same amplitude. However, our dark matter results yield a smaller bias 
parameter compatible with b ~ 2 on a scale of 20Mpc in contrast to previ- 
ous work which yielded to large bias factors. Our conclusions are thus more 
positive. According to the aspects analyzed in this work, global topological 
defect induced fluctuations yield viable scenarios of structure formation and 
do better than standard CDM on large scales. 



PACS numbers: 98.80-k 98.80.Hw 98.80C 



1 Introduction 



The formation of cosmological structure in the universe, inhomogeneities in the 
matter distribution like quasars at redshifts up to z ~ 5, galaxies, clusters, super 
clusters, voids and walls, is an outstanding basically unsolved problem within the 
standard model of cosmology. 

At first sight it seems obvious that small density enhancements can grow suf- 
ficiently rapidly by gravitational instability. But global expansion of the universe 
and radiation pressure counteract gravity, so that, e.g., in the case of a radiation 
dominated, expanding universe no density inhomogeneities can grow faster than log- 
arithmically. Even in a universe dominated by pressure- less matter, cosmic dust, the 
growth of density perturbations is strongly reduced by the expansion of the universe. 

On the other hand, we know that the universe was extremely homogeneous and 
isotropic at early times. This follows from the isotropy of the 3K Cosmic Microwave 
Background (CMB), which represents a relic of the plasma of baryons, electrons 
and radiation at times before protons and electrons combined to hydrogen. After a 
long series of upper bounds, measurements with the COsmic Background Explorer 
satellite (COBE) have finally established anisotropies in this radiation at the level 
of 

AT 

(— (0)) ~ KT 5 on angular scales 7° < 9 < 90° . 

On smaller angular scales the observational situation is at present somewhat con- 
fusing and contradictory [[|, but many upper limits require AT/T < 4 x 10~ 5 on all 
scales 9 < 8°. 

All observations together clearly rule out the simplest model of a purely baryonic 
universe with density parameter Q ~ 0.1 and adiabatic initial fluctuations (either 
the initial perturbations are too large to satisfy CMB limits, or they are too small 
to develop into the observed large scale structure). 

The most conservative way out, where one just allows for non-adiabatic initial 
perturbations (minimal isocurvature model), also faces severe difficulties || £|, ||, 
|6| . In other models one assumes that initial fluctuations are created during an 
inflationary epoch, but that the matter content of the universe is dominated by hot 
or cold dark matter or a mixture of both. Dark matter particles do not interact 
with photons other than gravitationally and thus induce perturbations in the CMB 
only via gravitation. In these models, inflation generically leads to Q — 1, while the 
baryonic density parameter is only f^/i 2 ~ 0.02, compatible with nucleosynthesis 
constraints. With one component of dark matter, these models do not seem to agree 
with observations 0, [7| , however, if a suitable mixture of hot and cold dark matter 
is adopted, the results from numerical simulations look quite promising p|, [TO 



although they might have difficulties to account for the existence of clusters at a 
redshift z ~ 1 |p"T . 



In these dark matter models initial fluctuations are generated during an inflation- 
ary phase. Since all worked out models of inflation face difficulties (all of them have 
to invoke fine tuning to obtain the correct amplitude of density inhomogeneities), 
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we consider it very important to investigate yet another possibility: Density pertur- 
bations in the dark matter and baryons might have been triggered by seeds. Seeds 
are an inhomogeneously distributed form of energy which makes up only a small 
fraction of the total energy density of the universe. Particularly natural seeds are 
topological defects. They can form during symmetry breaking phase transitions in 
the early universe [12|, |IBfl . Depending on the symmetry being gauged or global, the 
corresponding defects are called local or global. 

The fluctuation spectrum on large scales observed by COBE is not very far from 
scale invariant [Rj] . This has been considered a great success for inflationary models 
which generically predict a scale invariant fluctuation spectrum. However, as we shall 
see, also models in which perturbations are seeded by global topological defects yield 
scale invariant spectra of CMB fluctuations. To be specific, we shall mainly consider 
texture, 7r 3 -defects which lead to event singularities in four dimensional spacetime 
|T5| , |TB|] . Global defects are viable candidates for structure formation, since the 



scalar field energy density, ps, of global topological defects scales like ps oc l/(at) 2 
(up to a logarithmic correction for global strings) and thus always represents the 
same fraction of the total energy density of the universe (t is conformal time). 

ps/p ~ 8irG V 2 = 2e , (1) 

where r\ determines the symmetry breaking scale (see Fig. 1). For the background 
spacetime we assume a Friedmann-Lemaitre universe with = 1 dominated by cold 
dark matter (CDM). We choose conformal coordinates such that 

ds 2 = a 2 (—dt 2 + 5ijdx l dx^) . 

Numerical analysis of CMB fluctuations from topological defects on large scales has 
been performed in |17|, 18); a spherically symmetric approximation is discussed in 



||19|| . Results for intermediate scales angular are presented in [pOf] . All these inves- 
tigations (except [|TUJ) use linear cosmological perturbation theory in synchronous 
gauge and (except |18[) take into account only scalar perturbations. Here we derive 
a fully gauge invariant and local system of perturbation equations. The (non-local) 
split into scalar, vector and tensor modes on hyper surfaces of constant time is not 
performed. We solve the equations numerically in a cold dark matter (CDM) uni- 
verse with global texture. In this paper, we detail the results outlined in a previous 



letter [21]. Furthermore, we present explicit derivations of the equations, a descrip- 
tion of our numerical methods and we briefly discuss some tests of our codes. Since 
there are no spurious gauge modes in our initial conditions, there is no danger that 
these may grow in time and some of the difficulties to choose correct initial con- 
ditions (see e.g. fl8fl ) are removed. However, as we shall discuss in Section 3, the 
results do depend very sensitively on the choice of initial conditions. 

Nevertheless, we should keep in mind, that we are investigating models of struc- 
ture formation which rely on the particle physics and cosmology at temperatures 
of T ~ Tqut ~ 10 w GeV. An energy scale about which we have no experimental 
evidence whatsoever. The physical model adopted for our calculations should thus 
always be considered as a toy model, of which we hope it captures the features rel- 
evant for structure formation of the 'realistic physics' at these energies. Therefore, 
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we suggest, not to take the results serious much beyond about a factor of two or 
so. On the other hand, our models show that the particle physics at GUT scale 
may have left its traces in the distribution of matter and radiation in the present 
universe, yielding the exciting possibility to learn about the physics at the highest 
energies, smallest scales, by probing the largest structures of the universe. 

We calculate the CMB anisotropics on angular scales which are larger than the 
angle subtended by the horizon scale at decoupling of matter and radiation, 9 > 9&. 
For Q = 1 and z d « 1000 

9 d = « 0.03 « 2° . (2) 

It is therefore sufficient to study the generation and evolution of microwave back- 
ground fluctuations after recombination. During this period, photons stream freely, 
influenced only by cosmic gravitational redshift and by perturbations in the gravi- 
tational field (if the medium is not re- ionized). 

In Section 2 we derive a local and gauge invariant perturbation equation to calcu- 
late the CMB fluctuations. In Section 3, we put together the full system of equations 
which has to be solved to investigate gravitationally induced CMB fluctuations and 
the dark matter perturbation spectrum in a model with global topological defects. 
We discuss the choice of initial conditions and the numerical treatment of this sys- 
tem in Section 4. The next section is devoted to the presentation and analysis of 
our numerical results. We end with conclusions in Section 6. 

Notation: We denote conformal time by t. Greek indices run from to 3, 
Latin indices run from 1 to 3. The metric signature is chosen ( — h ++). We set 

k = C= k B oltzmann = 1 throughout. 



2 A Local and Gauge Invariant Treatment of the 
Perturbed Liouville Equation 

Collision-less particles are described by their one particle distribution function which 
lives on the seven dimensional phase space 

V m = {(x,p) E TM\g{x)(p,p) = -m 2 } . 

Here M. denotes the spacetime manifold and TAi its tangent space. The fact that 
collision-less particles move on geodesies translates to the Liouville equation for the 
one particle distribution function, /. The Liouville equation reads P2 



Xg(f) = . (3) 

In a tetrad basis (e^)^ =0 of Ai, the vector field X g on V m is given by (see e.g. p2[ ) 
X 9 = (p% - ^»P^) , (4) 
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where u v are the connection 1-forms of (Ai,g) in the basis e M , and we have chosen 
the basis 

d 

(e M )J =0 and (— )f =1 on TV m , p = p% . 

We now apply this general framework to the case of a perturbed Friedmann universe. 
The metric of a perturbed Friedmann universe with density parameter Q — 1 is given 
by ds 2 = g iiv dx ii dx v with 

9iw = a 2 {v^ + VO = a2 9^ , (5) 

where (rj^) = diag(—, +,+,+) is the flat Minkowski metric and (h^) is a small per- 
turbation, \h^ u \ <C 1. We now use the fact that the motion of photons is conformally 
invariant: 

We show that for massless particles and conformally related metrics, 

_ 2 ~ 

g^iv a g^-v i 

(X g f)(x,p) =0 is equivalent to (X§f)(x,ap) = . (6) 

This is easily seen if we write X g in a coordinate basis: 

8 

X — u*Fi — v 1 yp-yP— 



with 



— ~9 l ^(.gan,(3 +g/3^,a -gap,^ ) ■ 



2" 

The 6 M are the components of the momentum p with respect to the coordinate basis: 
P = P >±e ^ = b^df, . 

If (e M ) is a tetrad with respect to g, then e M = ae M is a tetrad basis for g. Therefore, 
the coordinates of of ap = ap^e^ = a 2 p fJ, e fJi = a 2 b^d^ with respect to <9 M on (M,g) 
are given by a 2 b^. In the coordinate basis thus our statement Eq. (|6]) follows, if we 
can show that 

(X g f)(x»,a 2 V) = iff (X g f)(x»,b i ) = (7) 

Setting v = ap = v M e M = w^d^, we have v M = ap M and w^ 1 = a 2 b^. Using p 2 = 0, we 
obtain the following relation for the Christoffel symbols of g and g: 

r a pb a b p = r a/3 b a b p + ^b a b i . 

For this step it is crucial that the particles are massless! For massive particles the 
statement is of course not true. Inserting this result into the Liouville equation we 
find 

a*X g f = W »(d,f\ b - 2^g) - f >V , (8) 

4 



where d^flb denotes the derivative of / w.r.t. x M at constant (b l ). Using 
d»f\ b = d,f\ w + 2^V^- : 

we see, that the braces in Eq. (R) just correspond to dfj,f\ w . Therefore, 

df 

a 2 X g f(x,p) = w^f\ w - r af3 w a w^^- = X g f(x,ap) . 

We have thus shown that the Liouville equation in a perturbed Friedmann uni- 
verse is equivalent to the Liouville equation in perturbed Minkowski space, 

(X s f)(x,v) = 0, (9) 

with v = f M e M = ap M e M .[] 

We now want to derive a perturbation equation for Eq. (§J). If e M is a tetrad in 
Minkowski space, e M = e M + \h v e v is a tetrad w.r.t the perturbed geometry g. For 
(x, v M e M ) G Pq, thus, (x, u M e M ) G Pq. Here Pq denotes the zero mass one particle phase 
space in Minkowski space and Pq is the phase space with respect to g, perturbed 
Minkowski space. We define the perturbation of the distribution function F by 

f{x, u%) = f(x, u%) + F{x, u%) . (10) 

Liouville's equation for / then leads to a perturbation equation for F. We choose 
the natural tetrad 

H = 9„ - \k Q » 
with the corresponding basis of 1-forms 

& = dx» + -h$dx u . 
Inserting this into the first structure equation, dQ^ = —u)^ v A dx u , one finds 

Using the background Liouville equation, namely that / is only a function of v = ap, 
we obtain the perturbation equation 

(d t + j l d l )F = ~[(h i0 - h mn )Y + (k, - h oj ,i )t¥]£ , 



■"■Note that also Friedmann universes with non vanishing spatial curvature, K ^ 0, are con- 
formally flat and thus this procedure can also be applied for K ^ 0. Of course, in this case the 
conformal factor a 2 is no longer just the scale factor but depends on position. A coordinate trans- 
formation which transforms the metric of K ^ Friedmann universes into a conformally flat form 
can be found, e.g., in Q. 
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where we have set v % = vy\ with v 2 = X^ =1 (V) 2 . Let us parameterize the perturba- 
tions of the metric by 

(/v)= (# t 2H L 6 tJ + 2H tJ )' 
with HI = 0. Inserting this above we obtain 

(d t + Yd t )F = -[H L + (A,i +\b i )! 1 + (Hij - Is y ) 7 V]t;^ . (12) 



From Eq. (|T2| ) we see that the perturbation in the distribution function in each spec- 
tral band is proportional to vjf. This shows ones more that gravity is achromatic. 
We thus do not loose any information if we integrate this equation over photon 
energies. We define 

7T r o , 

m = j / Fv dv . 

p r ar J 

4m is the fractional perturbation of the brightness l, 

i = a™ 4 J fv 3 dv . 
This is obtained using the relation 

4tt / %v A dv = -4 / fv 3 dvdtt = -Ap r a A . (13) 
J dv J 

Setting i = t(T(7, x)), one finds that i = (7r/60)T 4 (7, x). Hence, m corresponds to 
the fractional perturbation in the temperature, 

T( 7 ,x) =f (l + m( 7 ,x)) . (14) 

Another derivation of Eq. (|T4]) is given in ||24j| . Since the v dependence of F is of 
the form vj£, we have with Eq. ( |i~3|) 

F(x M , 7*, f ) = -m(x^,7>^ . 

This shows that m is indeed the quantity which is measured in a CMB anisotropy 
experiment, where the spectral information is used to verify that the spectrum of per- 
turbations is the derivative of a blackbody spectrum. Of course, in a real experiment 
located at a fixed position in the Universe, the monopole and dipole contributions to 
m cannot be measured. They cannot be distinguished from a background component 
and from a dipole due to our peculiar motion w.r.t. the CMB radiation. 

Multiplying Eq. (|T2|) with v 3 and integrating over v, we obtain the equation of 
motion for m 

d t m + 7 ^m = H L + (A n +~5 i ) 7 ' + (Hij ~ \b uj ) 7 V • (15) 
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It is well known that the equation of motion for photons only couples to the 
Weyl part of the curvature (null geodesies are conformally invariant). The r.h.s. of 
Eq. (|T5|) is given by first derivatives of the metric only which could at most represent 
integrals of the Weyl tensor. To obtain a local, non integral equation, we thus rewrite 
Eq. (|1|) in terms of Am. It turns out, that the most suitable variable is however 
not Am but \i which is given by 

X = Am - (AH L - |j=T,g ) - ^AB t - 3d^ tj )Y , 

where = -~{B hj +B jri ) + ^5i 3 Bf + i% 

Note that \ and Am only differ by the monopole contribution, AHl — (l/2)iJ* J ,j J ' 
and the dipole contribution, (1/2) (ABi — 3<9 J o"jj)7\ The higher multipoles of \ and 
Am agree. An observer at fixed position and time cannot distinguish a monopole 
contribution from an isotropic background and a dipole contribution from a peculiar 
motion. Only the higher multipoles, / > 2 contain information about temperature 
anisotropics. For a fixed observer therefore, we can identify A~ l \ with 5T/T. 

In terms of metric perturbations, the electric and magnetic part of the Weyl 
tensor are given by (see, e.g. [^5|, |24j) 

Eij = ^[AyiA-Hj-ay-AHij 

B{j ~ \^ilm^jmjl ~^~£jlm@imil ) > 

with Aij = didj — (l/3)5jjA . 

Explicitly working out (d t + Ydi)x using Eq. (|T5|) , yields after some algebra the 
equation of motion for \: 

(d t + J%)X = 3fa% + 7 VeM^% = S T (t, x, 7) , (18) 

where e^u is the totally antisymmetric tensor in three dimensions with £123 = 1. 
The spatial indices in this equation are raised and lowered with 5^- and thus in- 
dex positions are irrelevant. Double indices are summed over irrespective of their 
positions. 

In eqn. Eq. ( |18D the contribution from the electric part of the Weyl tensor does 
not contain tensor perturbations. On the other hand, scalar perturbations do not 
induce a magnetic gravitational field. The second contribution to the source term in 
Eq. ( |18D thus represents a combination of vector and tensor perturbations. If vector 
perturbations are negligible, the two terms on the r.h.s of Eq. ( |TBD yield a split into 
scalar and tensor perturbations which is local. 

Since the Weyl tensor of Friedmann Lemaitre universes vanishes, the r.h.s. of 
Eq. (pl|) is manifestly gauge invariant (this is the so called Stewart-Walker lemma 
|27j). Hence also the variable x is gauge invariant. Another proof of the gauge 



- 2 -Hi™5 l3 + Hi 3 +Hf l , i ] (16) 

(17) 
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invariance of Xi discussing the behavior of F under infinitesimal coordinate trans- 
formations is presented in |24[ . 

The general solution to Eq. ( |l8f) is given by 

X (t, x, 7) = f S T {t', x+{t'- t)j, 7 )df + xiU, x + (U- t)y, 7) , (19) 
Jti 



where St is the source term on the r.h.s. of Eq. Qlq). Let us compare this result with 
the more familiar one, where one calculates 5T/T by integrating photon geodesies 
(which is of course equivalent to solving the Liouville equation). For simplicity, 
we specialize to the case of pure scalar perturbations (the expressions for vector 
and tensor perturbations given in |24| can be compared with Eq. ( |I9"D in the same 
manner.) For scalar perturbations, integration of photon geodesies yields |24] 



(t f , x f , n) = + Vi ■ n* + (* - $)]{ + - $)dA . (20) 



ST 
~T 

Here \I/ and $ denote the Bardeen potentials as defined, e.g., in []26|, |24jl . On super 
horizon scales (which are the important scales for the Sachs-Wolfe) contribution 
Vi ■ n l can be neglected. Furthermore, the contributions in the square bracket of 
Eq. ( ]20|) from the final time t = tf, only lead to uninteresting monopole and dipole 
terms. We now use that the electric contribution to the Weyl tensor for purely scalar 
perturbations is given by |24|) 



E tJ = - iA)(* - *) = ±Ay(* - $) 

Therefore <9;(\l/ — $) = 38^ E^. Using furthermore 

|f = - £ f [i> - $ + (* _ , 

Eq. © leads to 

5 -^(t,x,n) = ±Dp(ti,Xi) -3j' A-WEijn'dt (21) 

If we take into account that the direction n in Eq. fl20|) , the direction of an incoming 
photon corresponds to —7 in Eq. (fT^), we find that Eq. ( PU[ ) coincides with Eq. ( TJ~9| ) 
for scalar perturbations, and that 

X(ti,^,7) = ^Mft,^) = ~^AD%\t i ,x- 1 {t-t i )) . (22) 

We now want to investigate this initial value and decompose Eq. ([H]) into terms due 
to CDM and terms coming from the source, the scalar field. We assume that dark 
matter and radiation perturbations are adiabatic on superhorizon scales, 



Dt ] = (4/3)£>J c > . 



S 



Since radiation and CDM probably have been a single fluid at early times (e.g. at 
the time of the phase transition), this assumption is reasonable. It is however, in- 
consistent to set = 4/3D( c > on subhorizon scales. Due to the different equations 
of state for the two components, adiabaticity cannot be maintained on sub-horizon 



scales |26[ . We can then derive from the equations given in |24 



\dP = h c + hc/(a/a) + $ 5 . 

Here the Bardeen potentials are split into parts due to cold dark matter (c) and the 
scalar field (5) respectively. For cold dark matter — —&c- Using this, we can 
bring Eq. ( |19|) into the form 

5T 1 2 • rf ■ 

— (t f , Xf ,n) = ~* c {t i ,x i )--Vc/{a/a){t i ,x i )+2J ^ c dt 

+& s (t l ,x l ) -3 J A^Srsfaxf - (t f -t)n,n)dt , (23) 
where Sts denote the portion of the source term due to the scalar field only: 

S TS = -M&Ef + n k nh klj dMp • (24) 



Eq. ( p3|) is much better suited for numerical investigation than the general expression 
Eq. ([19|). This can be demonstrated by considering the case of pure CDM without 
source term: In this case $c — — constant and from Eq. (|23| ) we easily recover 

the well-kown result 

5T 1 

— (t, x, n) = -^ c {ti, x-n(t- ti)) , 
whereas Eq. (|T9|) in this case leads to 



ST . 5T . . . 

— (t, x, n) = —{ti, Xi, n) + 2W £7 (t i , . 

In other words, the unknown initial condition in Eq. (|T9|) cancels 5/6 of the naive 
result for the case of adiabatic CDM fluctuations. Even though due to the existence 
of terms, the cancelation is slightly less substantial in our case, the assumption 
of adiabaticity on superhorizon scales is a crucial ingredient of the model. 

The electric and magnetic parts of the Weyl tensor are determined by the per- 
turbations in the energy momentum tensor via Einstein's equations. We assume 
that the source for the geometric perturbations is given by the scalar field and dark 
matter. The contributions from radiation may be neglected. Furthermore, vector 
perturbations of dark matter (which decay quickly) are neglected. The divergence 
of Ey is then determined by (see Appendix A) 

WE l3 = %nG Pc a 2 D l + 8irG(di6T 00 + 3(-)5T 0l - (3/2)<9%) , (25) 
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where the first term on the r.h.s. is the dark matter source term, pc denoting the 
dark matter energy density. The second contribution is due to the scalar field: The 
energy momentum tensor of the scalar field 



1 K 



A 



yields 

n . ^ T .. _ (a 2 /3)^-T/ = = <p n 4> d -(1/3)%(V0) 2 , 

ST oj = T 0j = M>,j , 

^oo = T o 5 o = i((0) 2 + (V0) 2 ), 

and Dj is a gauge invariant perturbation variable for the density gradient. For 
scalar perturbations Dj = djD. The evolution equation for the dark matter density 
perturbation is given by (see 



D + (~)D - 4nGa 2 p c D = 8ttG0 2 . (26) 
a 

During the radiation dominated era QuGprDr in principle has to be included in 
Eq. fl2"6j). But since radiation perturbations quickly decay on sub- horizon scales, and 



since dark matter fluctuations cannot grow in a radiation dominated universe f28 
their influence is not relevant. 

The equation of motion for is more involved. A somewhat cumbersome 
derivation (see Appendix A) yields 

a^iaBij)- - ABij = %itGS\f ] , (27) 

With Sif = —£lm(i8T t,j) m +£t m (ifj)i, m . 

Here denotes symmetrization in the indices i and j. 

To the source term only vector and tensor perturbations contribute. It is 
thus entirely determined by the energy momentum tensor of the scalar field. 



3 The System of equations for Global Scalar Field 
Induced Fluctuations 

In this section we collect all the equations which determine the system under con- 
sideration. We also repeat equations which have already been derived in Section 2. 
Let us begin with the scalar field equation of motion. 

The energy momentum tensor of the scalar field is a small perturbation. In 
first order perturbation theory, we can thus solve the equation of motion of the 
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scalar field in the background, Friedmann-Lemaitre geometry, neglecting geometric 
perturbations. The equation of motion for the scalar field (f) is then given by 

dV 

g^V,V u cj ) + — = 0, (28) 

where g^ v denotes the unperturbed metric and V M is the covariant derivative with 
respect to this metric. For our numerical computations, we consider an 0(4) model. 
In 0(N) models the scalar field, <fi G R N and the zero temperature potential is given 
by Vq = \{4> 2 — if) 2 for some energy scale r\. At high temperatures, T > T c ~ i], one 
loop corrections to the effective potential dominate and the minimum of the effective 
potential is at = 0. Below the critical temperature the minimum is shifted (in the 
simplest case) to < (ft 2 >= (1 — (T/T c ) 2 )t] 2 (see |12|, |16|] and references therein). The 



vacuum manifold, i.e. the space of minima of the effective potential, then becomes 
a (N — l)-sphere, 5 (Ar_1) . Since 



k < m 



the lowest non- vanishing homotopy group of a m-sphere is always ir m . Since proba- 
bly higher defects are unstable and decay into lower onesf], the m-sphere is a suitable 
vacuum manifold to study 7T TO defects. 

If the system under consideration is at a temperature T much below the critical 
temperature, T -C T c , it becomes more and more improbable for the field <ft to leave 
the vacuum manifold. <fi will leave the vacuum manifold only if it would otherwise 
be forced to gradients of order (V0) 2 ~ A0 2 r^ 2 , thus only over length scales of order 
/ = 1/(V%?) = m^ 1 (/ is the transversal extension of the defects). For GUT scale 
phase transitions / ~ 10 _30 cm as compared to cosmic distances of the order of Mpc 
~ 10 24 cm. If we are willing to loose the information of the precise field configuration 
over these tiny regions, it seems well justified to fix to the vacuum manifold M. 
Instead of discussing the field equation Eq. (pHF), we require G S {N ~ l) . The 
remaining field equation, □</> = 0, then demands that 

<j)/r) = (3 : M -> S {N ~ l) 

is a harmonic map from spacetime Ai into S^ -1 ^. 

The topological defects we are interested in are singularities of these maps. When 
the gradients of <fi become very large, like, e.g., towards the center of a global 
monopole, the field leaves the vacuum manifold and assumes non vanishing potential 
energy. If (3 G 5 (7V " 1} is enforced, a singularity develops by topological reasons. 

In the physics literature harmonic maps are known as a-models. The action of 
a cr-model is given by 

S a = [ 9^d^ A d v (3 B lAB ((3)J\g~\d 4 x , (29) 



2 This is an unproven conjecture, motivated, e.g., by observations of the density of textures and 
monopoles in liquid crystals and by numerical experiments [jl3| Q 
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where ^ab denotes the metric on S N ~ X and g^ u is the metric of spacetime. We now 
fix (3 to lay in the vacuum manifold, S N ~ X by introducing a Lagrange multiplier. 
We then obtain the following equation of motion for (3: 

□/?-(/?• □/?)/? = 0, (30) 

which shows that the cx-model is scale free. There are thus two possible evolution 
equations for the scalar field at low temperature. We call Eq. the 'potential 
model' evolution equation and Eq. (^) the a-model approach. 

The energy momentum tensor of the scalar field perturbs spacetime geometry 
and induces perturbations in the dark matter energy density according to Eq. ( f2~6l) 

D + {-)D - AixGa 2 p c D = 8ttG0 2 , (31) 
a 

where D is a gauge invariant variable for the dark matter perturbations |j28| . On 
subhorizon scales D ~ Sp/p. In comoving coordinates, the total perturbed energy 
momentum tensor is given by 

As already mentioned in section 2, the perturbed Einstein equations to this energy 
momentum tensor yield an algebraic equation for the divergence of the electric part 
of the Weyl tensor and an evolution equation for the magnetic part of the Weyl 
tensor (see Appendix A): 

$E l3 = -^Gp c a 2 A - 87tG(^W o + (-)5T 0l + ±3%-) , and (32) 
-{aB^y + -AB i3 = 8ttG4 S) , (33) 



with S\f ] = ej m(i [T f, i)m +f ,,■);, m ] , and r i3 - = <f),i<P,j ~$ij{^4>) 2 



The source term for the perturbation of Liouville's equation is given by Eq. (E 



- W&Ef + nVe^d^f = S ST {t, x, n) . (34) 

The CMB fluctuations are then determined according to 
ST 



T 



(t,x,n)= I Ssrit' ,x + {t' -t)n,n)dt' + ® s (ti, x + (t; - t)n, n) 
Ju 

+l^c(ti, Xi) - ^ifa/(a/a)(ti, a*) + 2 f ^ c dt (35) 



Eqs. (pS|) and ( pl|) to ( pop form a closed, hyperbolic system of partial differential 
equations. Actually all except the scalar field equation Eq. ( P8f) are linear pertur- 
bation equations with source term and can thus be solved, e.g., by the Wronskian 
method, i.e., by some integrals over the source term. The corresponding solution for 
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ST/T is given above in Eq. fl35|), the general solution of the dark matter equation is 
given below in Eqns. ((H), ([|6|) and (|47]) . 

Let us briefly describe the general solution for Bif We switch to Fourier space, 
because there the A is a simple multiplication by —k 2 and Eq. ( |3"BD becomes an 
ordinary differential equation with scalar homogeneous solutions 

b ± = -exp(±ikt) . (36) 

Ct 

The general solution to the inhomogeneous equation is given by 

B l] = (b + C± + b~C7 j ) + B^ m) , (37) 
where B hom denotes an arbitrary homogeneous solution and C + , C~ are given by 

r S (B) b- 

C± = -toGj-^-dt (38) 

C- 3 = SnGj^^dt. (39) 

Here W denotes the Wronskian determinant of the solutions which amounts to 

9ik 

W = b + b- - b-b + = q£ . (40) 

a 2 

4 Initial Conditions and Numerical Methods 
4.1 The scalar field: 

As already shown in the previous section, the equation of motion of the scalar field 
is given by 

^V M V^ + ^ = 0, (41) 
where g^ u is the background (unperturbed metric). With (3 = <f>/r] and m = y/Xr/, 



Eq. (5T) yields for O(N) models in a Friedmann universe 

d 2 t p + 2(a/a)d t p - V 2 (3 = ^a 2 m 2 {(3 2 - 1)0 . (42) 

This equation as it stands can not be treated numerically in the regime which is 
interesting for large scale structure formation. The two scales in the problem are 
the horizon scale t ~ (a/a) -1 and the inverse symmetry breaking scale, the comoving 
scale (am) -1 . At recombination, e.g., these scales differ by a factor of about 10 53 
and can thus not both be resolved in one computer code. 

There are two approximations to treat the scalar field numerically. As we shall 
see, they are complementary and thus the fact that both approximations agree with 
each other within about 10% is reassuring. The first possibility is to replace (am) -1 
by w, the smallest scale which can be resolved in a given simulation, typically twice 
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the grid spacing, w ~ 2Ax. The time dependence of (am) -1 which results in a 
steepening of the potential is mimicked by an additional damping term: 2(d/a) — > 
act/a, with a ~ 3 |JT] . Numerical tests have shown, that this procedure, which 



usually is implemented by a modified staggered leap frog scheme is not very 
sensitive on the values of a and w chosen. With this method we have replaced the 
growing comoving mass am by the largest mass which our code can resolve. For a 
(256) 3 grid which simulates the evolution of the scalar field until today, we obtain 
256Ax ~ to ~ 4x 10 17 sec/ao, so that w ~ 4x 10 15 sec/ao, i.e., am ~ r]a rec ~ 10 17 GeV 
is replaced by about w' 1 = a 10~ 39 GeV ~ 10~ 35 GeV, where we set a eq = 1. 

We are confident that this modified equation mimics the behavior of the field, 
since the actual mass of the scalar field is irrelevant as long as it is much larger 
than the typical kinetic and gradient energies associated with the field which are 
of the order the inverse horizon scale. Therefore, as soon as the horizon scale is 
substantially larger than Ax, the code should mimic the true field evolution on 
scales larger then w. But, to our knowledge, there exists no rigorous mathematical 
approximation scheme leading to the above treatment of the scalar field which would 
then also yield the optimal choice for a. 

Alternatively, we can treat the scalar field in the a-model approximation given 
in the previous section. This approach is opposite to the one outlined above in which 
the scalar field mass is much too small, since the cx-model corresponds to setting 
the scalar field mass infinity. 

The a-model equation of motion cannot be treated numerically with a leap frog 
scheme, since it involves non-linear time derivatives. In this second order 

accurate integration scheme has been developed by varying the discretized action 
with respect to the field |I8| . 

The two different approaches have been extensively tested by us and other work- 
ers in the field, and good agreement has been found on scales larger than about 



3-4 grid sizes |34|, We have compared our potential code with the exact 



spherically symmetric scaling solution |32[] and with our old spherically symmetric 



cr-model code ||19|| . Outside the unwinding events which extend over approximately 
3 grid sizes, the different approaches agree within about 5%. This is very encourag- 
ing, especially since the two treatments are complementary: In the a-model, we let 
the scalar field mass m go to infinity. In the potential approach, we replace m by 
~ 1/Ax ~ 200/to ~ 200a /10 10 ?/ ~ 10" 35 GeV. 

The integration of the scalar field equation is numerically the hardest part of the 
problem, since it involves the solution of a system of nonlinear partial differential 
equations. A good test of our numerical calculations, next to checking the scaling 
behavior of ps, is energy momentum conservation of the scalar field, T.\f^ u = 0. 
Energy momentum conservation in the potential model, about 15%, is slightly worse 
than in the a-model, where it is about 5% (see Fig. 2). Therefore, the final results 
presented here are all obtained with the cr-model approach. Our checks lead us to 
the conclusion, that we can calculate the scalar field energy momentum tensor, which 
then is the source of dark matter and CMB fluctuations to an accuracy of about 
10%. The problem of choosing the correct initial condition may induce another 
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(systematical) error in our calculations which we hope to remain below 20%. Other 
sources of error are negligible. 



4.2 Dark matter 

Once the scalar field /3(x,t) is known, the dark matter perturbations can easily 
be calculated by either using the Wronskian method (see below) or some standard 
ordinary differential equation solver. We have performed both methods and they 
agree very well. For later use, we briefly describe the Wronskian method. We 
normalize the scale factor by 

t 1 
a = —(1 H — t/r) , with 
r 4 



r = l/^/(4nG/3)p eq 



teq 



2(V2-1) ' 



Here t eq denotes the time of equal matter and radiation density, p ra d(teq) = Pc(teq) — 
(l/2)p(t eq ). We have normalized a such that a eq = a(t eq ) = 1. Transformed to the 
variable a, the dark matter equation Eq. ( p6|) then yields 

d 2 D ^ 2 + 3a dD 3 D — 2e$ 2 (^ a ) 2 — (1 + a)S/r 2 (43) 

da 2 2a(l + a) da 2a(l + a) dt 

S = 2e$ 2 and e = AtcGt] 2 . 



The homogeneous solutions to this linear differential equation are well known BB 



D 1 = l + ?o, (44) 

D 2 = (l + ^a)[ln(^n±l^_ 3v ^TT]. (45) 

The general solution to Eq. (|43D is given by 

D(t) = c^D^t) + c 2 (t)D 2 (t) (46) 

with 

ci = - J (SD 2 /W)dt , c 2 = J (SDjW)dt . (47) 

d(2o-l) 2a -1 



W = D X D 2 - D X D 2 



a\/ a + 1 ar 



is the Wronskian determinant of the homogeneous solutions. The integrals Eq. fl4"7|) 
have to be performed numerically with 5* = e(5 2 . When discussing the initial con- 
ditions for D in subsection 4.4, we shall present an analytic approximation for the 
source term S. 



15 



4.3 The CMB anisotropies 

The CMB anisotropies are given by 

5T A-l 

up to monopole and dipole contributions which we disregard. Here x is a solution 
Eq. (|T9| ) of Eq. (|TJ). The source term St is determined via Eq. (|25|) and Eq. (p7|) . 



However, using this straightforward approach results in a big waste of computer 
memory (which we cannot afford): We would be satisfied to calculate ST/T for about 
30 observers in each simulations, which means we need A _1 x only at 30 positions 
x. But since we have to perform an inverse Laplacian which is done by fast Fourier 
transforms, we have to calculate x on the whole grid, which consists of 192 3 ~ 8 • 10 6 
positions. In Addition, to calculate the spherical harmonic amplitudes of ST/T up to 
about I ~ 40 (angular resolution of about 4°), we need typically 5 x 10 4 directions n. 
The x variable alone (in double precision) would thus require 16 G-bytes of memory, 
an amount which is not available on most present day computers. The way out is to 
take the inverse Laplacian already in the equation of motion Eq. fllS]) . This results 
in 

(dt + Yd^Y = -3f A-^'iy -tV'^A- 1 ^,^) = A- 1 ^,^) -(48) 

Here the inverse Laplacian has to be performed for a vector field and a symmetric 
traceless tensor field, a total of 8 scalar variables which only depend on x and not on 
-y. For a 192 3 grid this requires nearly 1 Gbyte of memory, no problem for presently 
available machines. Eq. (H3) has the general solution (see Eq. (|23|)) 



5T Q rf , 

— (to,aj 0> 7) = +® (t i ,x i )-3j A S T g(t,x — (tf — t)n,ri)dt 

+U> c (t i} x i ) - ^ifc/(a/a)(jt i ,x i ) + 2 J [ i> c dt , (49) 



with Sts given in Eq. (|24|). The first term of Eq. (^) determines the initial condition 



of the CMB anisotropies due to the source term. In the numerical simulation we just 
set it -3n l /\-\d j E^^x^U. This assumes that the source term is approximately 
constant until tj and that magnetic contributions can be neglected. The resulting 
amplitude is not very sensitive to this assumption, but changing it can somewhat 
change in the spectral index. We have solved Eq. (^) numerically by just summing 
up the contributions from each time step for 27 observer positions x . The value 
of the source term at position x + (t — to) is determined by linear interpolation. 
The quantity d l E^ is determined by Eq. ( p5|) and its inverse Laplacian is calculated 
by fast Fourier transforms. To obtain A _1 fi^- from equation Eq. (^7|), we directly 
calculate A _1 Sb in fc-space, then solve the ordinary, linear differential equation for 
/\^ 1 B i j in fc-space by the Wronskian method. Since, as we shall argue later, all 
components Tf£> in average scale like A/y/t on super horizon scales, S^B) ~ At~ 3 / 2 
and therefore C ± oc t 3//2 on super horizon scales. Therefore, we can neglect the 
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contribution to C from the lower boundary in the integral. Furthermore, since 
the homogeneous solution B^ ™^ is decaying, we drop it entirely. This procedure 
corresponds to setting B(ti) = and calculating B(t) according to Eq (]37|). 



4.4 Initial conditions 

Initially, the field itself and/or the velocities <fi are laid down randomly on the grid 
points. The initial time, ti n is chosen to be the grid size, U n = Ax, so that the field 
at different grid points should not be correlated. The configuration is then evolved 
in time with one of the approximation schemes discussed above. 

Because our initial conditions for dark matter and photons very sensitively de- 
pend on the scaling behavior of the dark matter, we can only start the dark matter 
or photon simulations when scaling is fully reached, ti n = 8Ax. Starting our simula- 
tions, e.g., at t — 4Ax changes the results by about a factor of 2. Further doubling 
of the initial time, changes our results by less than 20%, we thus believe that at 
t = 8Ax scaling is sufficiently accurately. Unfortunately, this late initial time re- 
duces our dynamical range to about 192/8 = 24 for a 192 3 grid, which is seen clearly 
in our results for the CMB anisotropics discussed below. 

It is very important to choose the correct initial conditions for the dark matter 
and, especially, photon perturbations. Changing them can change the CMB fluctua- 
tion amplitudes by nearly a factor two. Since these fluctuations are used to normalize 
the model, i.e. to determine e, this reflects in corresponding changes in e. We want 
to do better than a factor two by choosing physically plausible initial conditions. 
The cleanest way would by to simulate the evolution of perturbations through the 
phase transition, assuming that before the phase transition, the universe was an 
unperturbed Friedmann universe with = 0. On the other hand, since we want 
to calculate the perturbation spectrum on scales of up to lOOOMpc with a (256) 3 
grid, we cannot start our dark matter and CMB simulation earlier than at a time 
when the horizon distance is approximately 8Ax ~ 30Mpc. At the beginning of the 
scalar field simulation our grid scale Ax ~ 4Mpc is of the order of the horizon scale. 
We therefore have to decide on the amplitudes of super horizon perturbations. One 
possibility is setting all geometrical perturbations initially to zero. The requirements 

3%-(*i) = and Sfj(U)=0 

then yield initial conditions for the dark matter fluctuations D and the photon 
variable \. But these, let us call them 'strict isocurvature' initial conditions are not 
natural since they do not propagate in time: Even if we start with E and B vanishing 
on super horizon scales, after some time residual fluctuations have leaked into these 
scales and one obtains the white noise fluctuations spectrum on super horizon scales 
shown in Fig. 3. This does not violate causality, since white noise is uncorrelated 
and just results from the residuals of correlated fluctuations on smaller scales. The 
correct initial values for D and D would of course be those obtained by solving the 
equation of motion Eq. fl2"E|) from the symmetry breaking time until the start of the 
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simulation. We found a method to incorporate this at least approximately: The 
spectrum of the dark matter source term 87rG|0 2 | 2 can be approximated by 

^ = 2 ^ = 'JvJ***'** 53 WT^T^W) • (50) 

with 

A = 3.3 , ai = -0.7/(2tt) , a 2 = 0.7/(2vr) 2 . 

This numbers have been obtained by a ^-minimization scheme. The approximation 
is not very good. It yields a \ 2 ~ 2000 for about 1000 data points. Its comparison 
with the real data in Figs. 4 and 5 shows that Eq. ( |50|) approximates the source 
term to about 10% on superhorizon scales, but does not follow the wiggles present 
in the data on smaller scales. Since we shall not use the fit on subhorizon scales, 

this is not important for our simulations. However, in general <fi 2 is complex and 
setting it equal to its absolute value, we neglect the evolution of phases. Again, by 
causality, this will not severely affect scales larger than the horizon, since on these 
scales the phases are (approximately) frozen. But on subhorizon scales our fis is nor 
very useful due to the incoherent evolution of phases. Assuming this form of the 
source term, we can solve Eq. fl26|) analytically on super horizon scales, where we 
approximate the source term by 

2ej3 2 = -y=, on super horizon scales. (51) 



The homogeneous solutions of Eq. (^6|) are given by Eq. (|44| , fi5)) The general inho- 
mogeneous solution, D = C\D\ + c 2 D 2 , even with the simple source term Eq. (|5~TD , 
becomes rather complicated. But in the radiation and matter dominated regimes 
we find the simple approximations 

D = (A/7)t 2 S ; D = (<o/l)t 2 S radiation dominated (52) 
D = -(A/9)t 2 S ; D = -(2/3)t 2 S matter dominated. (53) 

From D we can calculate leading to the dark matter contribution to the 
CMB anisotropies. 

As mentioned above, the initial contribution of the scalar field is approximated 

by 

$ 5 (t i ,a; J )~-3^A- 1 (^.f ) . 

The result is very sensitive to initial conditions: If we do not separate the dark 
matter and choose some arbitrary, non-adiabatic initial condition, the resulting C/s 
increase by nearly a factor of 10 and the dark matter induces 80% of the total 
fluctuation. However, choosing the adiabatic initial condition discussed in Section 3, 
leading to Eq. (f^Up, dark matter only contributes about 20% to the C/s and the 
main contribution is due to the defects. The dark matter contribution to the CMB 
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anisotropies is not scale invariant, but white noise. It has spectral index n = 0. This 
result was found numerically (see Fig. 6) but, as we argue in Appendix B, it can 
also be understood analytically. 

our value of e obtained with these physical isocurvature and on super horizon 
scales adiabatic initial conditions is in reasonable agreement with the values obtained 
3 and |T7| . 



in 



Let us also present a heuristic derivation of the numerical finding Eq. (|5ll) on 
superhorizon scales: We know that the average value {(3 2 ) oc 1/t 2 , the usual scaling 
behavior. The Fourier transform of f3 2 determines the fluctuations on this 'back- 
ground' on a given comoving scale A = 2ii/k. As long as this scale is super horizon, 
A > t, a patch of size A 3 consists of N = (A/t) 3 independent horizon size volumes. 
The fluctuations on this scale should thus be proportional to 



(3 2 oc (/3 A )/VN oc 1/y/t , 

which is just the behavior which we have found numerically on super-horizon scales. 

As soon as a given scale becomes sub-horizon, A <C t, (5 2 starts decaying from 
this large scale value like 1/t 2 . 



5 Results 

5.1 CMB anisotropies 

To analyze the CMB anisotropies, we expand ST/T in spherical harmonics 
ST 

— (to, a, 7) = 53 ai m {x)Y lm (j) . (54) 

lm 

As usual, we assume that the average over N x different observer positions coincides 
with the ensemble average and define 

^C^EW-)!'. t>2. (55) 

Gaussian fluctuations are characterized by the two point correlation function. Since 
the angular two point correlation function is given by 

(^(n)^(n')) (n • n^e) = ^(2i + l)C e P e (cos6), (56) 

Gaussian distributed CMB fluctuations are fully determined by the Q's. However, 
as can be seen from Fig. 7, in our case the distribution of the CMB fluctuations is 
not quite Gaussian. It is slightly negatively skewed. We find an average skewness of 
—0.5 and a kurtosis of 0.7. For a simulation on a 192 3 grid with 27 different observer 
positions for each simulation. The harmonic amplitudes with are shown in Fig. 8. 
The low order multipoles depend strongly on the random initial conditions (cosmic 
variance), like in the spherically symmetric simulation [HIT]. 
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It is well known, that cold dark matter fluctuations with a power spectrum of 
spectral index n gravitationally induce CMB anisotropics with a spectrum given by 

m 

r(i + (n-l)/2)r((9-n)/2) 
1 2 r(/ + (5-n)/2)r((n + 3)/2) • 1 ) 

We have performed a least square fit of our numerical results for log(C|) and the 
log(CV) obtained from Eq. ( [STD -Fl If we take into account all the C/s reliably calcu- 
lated in our simulations, which limits us approximately to £ < 22 we find a very nice 
scale invariant spectrum, 

n = 0.9 ±0.2 (58) 

with quadrupole amplitude 



Q = ^{h/ATi)C 2 T CMB = 2.8 ± 0.7K ■ e . (59) 

The 1, 2 and 3 sigma contour plot is shown in Fig. 9. The minimal \ 2 is 0.56. 

It is very interesting, that the dark matter contribution to the CMB anisotropies 
does not yield a scale invariant spectrum, but white noise. This can be understood 
analytically: The contributions to 5T/T in Eq. (|23|) are not very important and 



, , m;mN , , , 1 T , , N . eA exr>(ik • nt ) 

(6T/T)c(t , k) ~ -* c (t f , k) exp(zfc • nt ) 



on super horizon scales. For the second equal sign we used k 2 ^> c = 4:7tGD c ~ 
By standard arguments (see, e.g. P7|| ) one then finds 

corresponding to Eq. (p^) with n = 0. This is also what we find numerically (see 
Fig. 6). The dark matter contribution caused the spectral index n of the total CMB 
anisotropies to drop slightly below n = 1. 

To reproduce the COBE amplitude Qcobe = (20 ± 5)fj,K [nj, we have to nor- 
malize the spectrum by choosing the phase transition scale r\ according to 

e = 47r(V= (0.8±0.4)10~ 5 . (61) 

This value is somewhat smaller, but still comparable with the results obtained in 



ni, ig. 

Another method to determine e is the following: The total temperature fluctua- 
tion amplitude on a given angular scale 9 C is given by 

4(0 C ) = E Ct(2£ + 1) exp(-£ 2 ^ 2 /2) . (62) 



In Fig. 10 we show ot as a function of 9 C . In a recent analysis of the COBE data[[39[ 

a^p OBE \7°) ~ AAfiK and a^ OB (10°) ~ 40fiK for a spectral index n ~ 1, which 

leads again to the result given in Eq. flST]) . 

3 In the case of topological defect induced fluctuations, the spectrum does not have exactly 
this form, since CBM fluctuations are not only induced by the dark matter but mainly by the scalar 
field perturbations and the assumptions made for the derivation of this formula are not valid. 
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5.2 Dark matter fluctuations 



Using fast Fourier transforms we calculate the spectrum P(k) = \S(k)\ 2 of the dark 
matter density fluctuations is shown in Fig. 11. The fit shown as solid line in Fig. 11 
is given by 

Ck 

= {l + ak + m , i+[lkfY , (63) 

with h = 0.5 and 

C = 215/i~ 1 Mpc 4 (64) 
a = l0h- 2 Mpc (=0.5r) (65) 
(3 = 1.25h~ 2 Mpc ~ r/(4vr) (66) 
7 = 2.3h- 2 Mpc ~ t/(2tt) , (67) 

where we have used r = 19.36/i~ 2 Mpc which is approximately the comoving time at 
equal matter and radiation. 

The parameter C, which is most important to determine the bias factor can 
also be obtained by the following rough analytical argument: On super horizon 
scales, \D\ 2 ~ (0.5eA) 2 t 3 according to Eq. (|53"D. As soon as the perturbation enters 
the horizon at t = 2n/k, the source term disappears and D starts growing like t 2 , 
leading to 

pfl , , , (0.5eA) 2 ; ^ 4 (2tt) 3 

P{k,t ) kt = ——C an k. (68) 

Inserting the numbers e = 0.8 x 10~ 5 , A = 3.3, = 4aor 2 , a ~ 2.5h 2 x 10 4 , we 
obtain C an ~ 190 in excellent agreement with Eq. (|64"D. Figure 11 can be compared 



directly with the IRAS observation ]3S] and it is compatible with a bias factor of 
order 1. A more detailed calculation with Gaussian or square hat window function 
yields for e = 0.8 x 10 -5 

a sim (10MPc) = I/610 -0.5-1 cx QDO t{10MPc) = 1 for h = 0.5, (69) 

yielding 610 ~ 1 — 2 for the value of e found by comparison with COBE Eq. (|6T| ) A 
value even somewhat closer to 1 is found for 620- Observations and simulations of 
nonlinear clustering of dark matter and baryons j|0 suggest a bias factor 6 10 ~ 1 — 2 



which is compatible with our results. It is remarkable, that unlike in the simula- 
tions by Pen et al. |18||, our bias factor is approximately constant and physically 



acceptable. (To determine our power spectrum, we have not taken into account any 
smoothing which might change the results by at most 15%.) In Fig. 12 we have 
shown the dark matter pixel distribution from a 100 3 simulation. It is interest- 
ing, that the skewness of the dark matter distribution is positive, where the ST/T 
skewness is negative. 
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6 Conclusions 



Our simulations show that global texture lead to a scale invariant spectrum of mi- 
crowave background fluctuations on large scales like inflationary models of structure 
formation. This is one of the main results of this investigation. It is however in- 
teresting that the dark matter contribution to the CMB anisotropics is not scale 
invariant, but white noise. Therefore it is important that the initial condition for 
dark matter and radiation are adiabatic in which case the dark matter contribution 
to the C/s is small and the flat spectrum caused by the defects is maintained. 

Our second main result is the dark matter fluctuation spectrum. The spectrum 
is very close to scale invariant and the bias factor needed for e from the CMB 
anisotropies is around b ~ 1-2. This value is certainly acceptable and smaller than 
the bias factor obtained in previous investigations [?]. 

The deviation from Gaussian statistics seems to us not very significant (see Figs. 7 
and 12) and it is thus important to develop other means to distinguish topological 
defects from inflationary scenarios. A clean and promising candidate for this dis- 



tinction are the Doppler peaks which are calculated for the texture scenario in Wl 

From our investigations we thus conclude that concerning the large scale CMB 
anisotropies and the linear dark matter perturbation spectrum the texture scenario 
and probably also other models with global defects are compatible with present 
observations. 
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A The equation of motion for the magnetic part 
of the Weyl tensor 

The Weyl tensor of a spacetime (Ji4,g) is defined by 

C^ p = W\ p - 2«] + , (Af ) 

where denotes anti-symmetrization in the indices /i and v. The Weyl curvature 

has the same symmetries as the Riemann curvature and it is traceless. In addition 
the Weyl tensor is invariant under conformal transformations: 

C%fo) = C^ p (a 2 g) 

(Careful: This equation only holds for the given index position.) In four dimensional 
spacetime, the Bianchi identities together with Einstein's equations yield equations 
of motion for the Weyl curvature. In four dimensions, the Bianchi identities, 

are equivalent to PB 



C apjS. g = jptefl _ _ g i[<* R A . (A2) 
6 

This together with Einstein's equations yields 

C a ^ s ; s = 87rG(T 7[a;/3] - -g^ a T^) , (A3) 

3 

where T pu is the energy momentum tensor, T = T^. 

Let us now choose some time-like unit vector field u, u 2 = —1. We then can 
decompose any tensor field into longitudinal and transverse components with respect 
to u. We define 

h\ = g\ + u»u v , 

the projection onto the subspace of tangent space normal to u. The decomposition 
of the Weyl tensor yields its electric and magnetic contributions: 

E pu = C^ va u x u° (A4) 

B,u = \c^u X r?ly- (A5) 

where r] al3lS denotes the totally antisymmetric 4 tensor with 770123 = \J—Q- Due 
to symmetry properties and the tracelessness of the Weyl curvature, E and B are 
symmetric and traceless, and they fully determine the Weyl curvature. One easily 
checks that E pu and B pu are also conformally invariant. We now want to perform 
the corresponding decomposition for the energy momentum tensor of the scalar field, 
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We define 



1 

3 



-TffhT (AT) 



In 



ps = 3j?ttV (A6) 
Ps = 

-h?T<gu* qi = --T f (A8) 

V = hfhj^g - hrfs ■ (A9) 

We then can write 

= PsUfiU u + p s h^ u + q^u u + u^q u + r Mi , . (A10) 

This is the most general decomposition of a symmetric second rank tensor. It is 
usually interpreted as the energy momentum tensor of an imperfect fluid. In the 
frame of an observer moving with four velocity u, ps is the energy density of the 
scalar field, ps is the isotropic pressure, q is the energy flux, u ■ q = 0, and r is the 
tensor of anisotropic stresses, r^ u h^ u = t^ v u^ = 0. 

We now want to focus on a perturbed Friedmann universe. We therefore consider 
a four velocity field u which deviates only in first order from the Hubble flow: u = 
(l/a)do+ first order. Friedmann universes are conformally flat, and we require the 
scalar field to be a small perturbation on a universe dominated by radiation and cold 
dark matter (CDM). The energy momentum tensor of the scalar field and the Weyl 
tensor are of thus of first order, and (up to first order) their decomposition does not 
depend on the choice of the first order contribution to u, they are gauge-invariant. 
But the decomposition of the dark matter depends on this choice. Cold dark matter 
is a pressure-less perfect fluid We can thus choose u to denote the energy flux of 
the dark matter, Tj^u v = —pcu^. Then the energy momentum tensor of the dark 
matter has the simple decomposition 

T^= Pc u,u v . (All) 

With this choice, the Einstein equations Eq. (|A3|) linearized about an Q = 1 Fried- 



(S) 

mann background with T^ kground = yield the following 'Maxwell equations' for E 
and 

i) Constraint equations 

d % B l3 = ^Grjjp^q^ (A12) 
tfEn = %^G^a 2 Pc D, j +^a 2 Ps , j -^r ij -^q j ) . (A13) 

ii) Evolution equations 

aBij + dBij-a^^^E^ 5 = -^Ga 2 h a{i r 1]) ^ u u P T a ^ (A14) 

Eij + ^Eij + ah^ a r] j ) f3l su /3 B a r ' s = -A-nG(aq i:j - -t^ + + apc%iAl5) 
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where denotes symmetrization in the indices i and j. The symmetric traceless 

tensor fields q^ v and u^ v are defined by 

1 



1 \ 

u fiu — u (tt;v) ~^ l nv u ; A • 

In Eqs. ( |A14jj and ( |A15| ) we have also used that for the dark matter perturbations 
only scalar perturbations are relevant, vector perturbations decay quickly Therefore 
u is a gradient field, Ui = U-i for some suitably chosen function U. Hence the vorticity 
of the vector field u vanishes, U[^. v ] = 0. With 

Voijk = a 4 e ijk , p s = a~ 2 T^ and = -a~ x T^ , 



we obtain from Eq. (| A 1 3|) 



d% = MG{ l -p c a 2 D )J +±T&,j ~d inj + -T^) . (A16) 



In Eq. ( |A16| ) and the following equations summation over double indices is under- 
stood, irrespective of their position. 

To obtain the equation of motion for the magnetic part of the Weyl curvature we 
take the time derivative of Eq. ( |A14| ), using u = (l/a)do + 1. order and rjoijk = a 4 %fe- 
This leads to 

(aBij)" = —a(ei m (i[Ej)i + -E j )i}, m -A7cGei m ( i [f j )i, m +-T j - ) i, m }) , (A17) 

where we have again used that u is a gradient field and thus terms like eijkUij,k vanish. 
We now insert Eq. (|A15|) into the first square bracket above and replace product 
expressions of the form ey^e^ m and e^e^^ with double and triple Kronecker deltas. 
Finally we replace divergences of B with the help of Eq. ( |A12|) . After some algebra, 
one obtains 

£im(i[Ej)i + —Efii], m = —AB^ — ^Gei m (i\2aqi m j) +fjy, m ~^2 T j)hm ] ■ 

Inserting this into Eq. (|A17 ) and using aqi = —T^ = — <f>(f>,i, we finally obtain the 
equation of motion for B: 

a _1 (a5)v. - AB^ = SttGS^ , (A18) 

with 

Sif } = eim(i[-Ta,j) m +Tj)i, m ] , and Uj = 0,; </>,j -^<%(V0) 2 . (A19) 

Since dark matter only induces scalar perturbations and B^ consists of vector and 
tensor perturbations, it is independent of the dark matter fluctuations. Equations 
Eqs. ( |A16| ) and (|A18 ) are used in section 2, where we need d i E i j and Bij as source 
terms in the Liouville equation. 
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Figure Captions 

Fig. 1 

The scaling behavior for (p+3p)a 2 found numerically in (128) 3 simulations for differ- 
ent O(N) models. The time is given in units of the grid spacing Ax. For comparison 
a dashed line oc 1/t 2 is shown. After some initial oscillations, for iV > 3 the scaling 
is very clean until; t ~ 80, where finite size effects can become important, clean 
until 

Fig. 2 

The spectrum of the electric part of the Weyl tensor as a function of kt at time t — 8 
for a grid of size 160. On large scales, kt/ (2ir) < 1, the spectrum is flat, white noise. 



Fig. 3 

The quantities |(T %)| 2 (dotted line), |(7f; M )| 2 (dashed lines) and |(T °/t)| 2 (solid 
line) are shown as functions of k. The inaccuracy in energy and momentum conser- 
vation is below 10% for k < 32 = 128/4. This lets us that our code is accurate to 
better than 10% for wavelengths of 4 grid spacings and larger. 

Fig. 4 

The dashed curves and the triangles show f3 2 as a function of t for fixed values of 
k = n/ritot for n tot = 128. The solid lines show the fits according to the fitting 
formula given in the text. 

Fig. 5 

The crosses and triangles show /3 2 as a function of n, with k = n/n tot for fixed values 
of t. The solid curves show the fits according to the fitting formula given in the text. 



Fig. 6 

The dark matter contribution to the Q's from a (160) 3 simulation. (£ + 2)(£ + 
l)£Ct/2A is shown. For I < 20 which is the dynamic range of this simulation, a 
white noise, n — 0, spectrum fits reasonably well. 

Fig. 7 

The pixel distribution of ST/T for one observer. 
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Fig. 8 

The values £{l + l)Q/6 for 27 observers are plotted for e = 1. The crosses are the 
individual observers and the solid line indicates the average. The sharp drop after 
£ ~ 30 is due to finite resolution (our dynamical range is approximately 25). 

Fig. 9 

The x 2 contour plots for 66%, 95% and 99% confidence levels from fitting the Ce, 
to a spectrum with index n with quadrupole amplitude Q according to Eq. ([)7]) for 
£ < 20. In total 81 observers from 3 different (192) 3 simulations have been taken 
into account. 

Fig. 10 

The root mean square of the temperature fluctuation at given angular scale is shown 
as a function of angle. 

Fig. 11 

The dark matter power spectrum (without bias and nonlinear evolution). The result 
is averaged over 15 simulations on a (256) 3 grid of two different physical scales. The 
error bar indicates one standard deviation. 

Fig. 12 

The dark matter pixel distribution from linear perturbation theory. The positive 
skewness (0.76) and positive kurtosis (1.2) are clearly visible. 
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